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ABSTRACT 

The Cygnus X-ray Emission Spectroscopic Survey (CyXESS) sounding rocket 
payload was launched from White Sands Missile Range on 2006 November 20 and 
obtained a high resolution spectrum of the Cygnus Loop supernova remnant in 
the soft X-rays. The novel X-ray spectrograph incorporated a wire-grid collima- 
tor feeding an array of gratings in the extreme off-plane mount which ultimately 
dispersed the spectrum onto Gaseous Electron Multiplier (GEM) detectors. This 
instrument recorded 65 seconds of usable data between 43-49.5 A in two promi- 
nent features. The first feature near 45 A is dominated by the He-like triplet of 
O VII in second order with contributions from Mg X and Si IX-Si XII in first 
order, while the second feature near 47.5 A is first order S IX and S X. Fits to 
the spectra give an equilibrium plasma at log(T) = 6.2 {kTe = 0.14 keV) and 
near cosmic abundances. This is consistent with previous observations, which 
demonstrated that the soft x-ray emission from the Cygnus Loop is dominated 
by interactions between the initial blast wave with the walls of a precursor formed 
cavity surrounding the Cygnus Loop and that this interaction can be described 
using equilibrium conditions. 

Subject headings: instrumentation: spectrographs — ISM: individual (Cygnus 
Loop) — line: identification — supernova remnants — X-rays: individual (Cygnus 
Loop) 



1. Introduction 



Supernovae greatly influence the dynamics within the interstellar medium (ISM). Their 
ubiquitous nature and the size of their remnants allow them to influence multiple phases of 
the ISM while influencing the evolution and structure of the galaxy. They are responsible 
for chemical enrichment of the ISM, play a major role in energy input, and are do minant 
sources for hot ionized gas, which fills much of the galaxy (IMcKee fc Ostriker Ill982l ). 
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One of the mo st studied SNRs is the Cygnus Loop. It is a nearby r emnant at 540 pc 



( iBlair et ah Ill975l ). quite large, filhng ~ 3° x 3° (jLevenson et al. 1119971 ). and is relatively 
unabsorbed {E{B - V) = 0.08; IPesen et al. 1 Jl982h l The blast wave is currently encoun- 
tering the surrounding ISM leading to bright emission at all wavelengths. Therefore, the 
Cygnus Loop serves as an excellent test bed for SNR evolution theory. A comprehensive 
and detailed imaging study of the remnant's X-ray and optical emission was performed by 
Levenson et al. I (119971 . Il998l ). One of the major findings of these studies is that the Cygnus 
Loop morphology is not indicative of the typical theoretical p icture of a blast wave propa- 
gating through a uniform medium (ISedov Ill993l : ISpitzer Ill998l ). but instead is an interaction 
of the blast wave with an inhomogeneous medium. The precursor formed cavity is nearly 
spherical and surrounded by high density clumps and lower density gas. This causes most 
emission to originate from a limb-brightened shell. The optical emission occurs when the 
shock is decelerated rapidly in high density clumps. The X-ray emission occurs in these 
clumps as the blast wave interacts with the cloud and hardens toward the interior as a re- 
flection shock propagates inward reschocking the already shocked gas. X-rays also originate 
as expected in the lower density gas that follows the theoretical model more closely. 



The Cygnus Loo p has been spectroscopically observed in the soft X-ray to some extent. 
Vedder et al. I (Il986l ) obtained a high resolution (E/AE~ 50) spectrum at energies below 1 
keV using the Focal Plane Crystal Spectrometer (FPCS) on the Einstein Observatory, which 
observed a 3' x 30' bright northern section midway between the center and the limb. They 
detected O VIII Lyman-a at 653 eV, the resonance line in the He-like triplet of O VII at 
574 eV and a blend of Ne IX lines between 905-920 eV. They do not detect the forbidden or 
intercombination lines of the O VII triplet. Even though the detected fluxes are consistent 
with an equilibrium plasma with solar abundances at T ~ 3 x 10^ K, the authors argue that 
the plasma is not yet in collisional ionization equilibrium (CIE) due to the lack of forbidden 
O VII emission. 



ROSA T Position Sensitive P r oport ional Counter (PSPC) observations by lLevenson et al 



(119991 ) an d iMivata fc Tsunemi I (120011) show results consistent with the previous imaging 
studies. A iLevenson et al. I (119991 ) spectrum of the decelerat ed cloud shock was best fi t by a 
Raymond- Smith equilibrium plasma at cosmic abundances (IRaymond fc Smith 1119771 ) with 
kTe = 0.06(+0.04, -0.02) keV and absorbing column Nh = 7 x W'^°{±4 x 10^°) cm'^. A 
spectrum extracted from a region interior to the cloud shock had a similar low temperature 
component, but also included a high temperature contribution, kTei = 0.25(+0.14, —0.06) 
keV, from the reflected shock, which hardens the interior emission. The low temperature 
corresponds to a shock velocity of Vg = 230 km s~^, which has slowed considerably from the 



initial blast wav e velocity of v,^ ~ 500 km s ^ (ILevenson et al. Ill998l ) resulting in an equi- 



librium plasma. iMiyata fc Tsunemi I (120011 ) also flnd that equilibrium plasmas with cosmic 
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abundances fit the low temperature components of tlieir extracted spectra {kT^ = 0.043 to 
0.067 keV) and tliat tliese contributions are greatest toward tlie limb. However, they use 
ASCA data and a nonequilibrium model with depleted abundances to fit the interior harder 
spectra [kT^ = 0.27 to 0.34 keV). 

The Chandra X-ray Observatory Advanced CCD Irnaging Spectrometer (ACTS) al so ob- 
served the Cygnus Loop (ILevenson et al. II2002I : iLeahy II2004I ). iLevenson et al. I (120021 ) used 



the Xspec ( Arnaud. K. A. 1996h spectra. 



MEKAL equilibrium model (IMewe et al. 



19851 : lArnaud fc Ravmond 1^99^ . 



fitting software package and more specifically the 



19851. 



Liedahl et al. 



986l: iKaastra 1119921 : lArnaud fc Rothenflug 



199 



to fit the extracted spectra. Again, 



the softer emission originates near the limb {kTe ~ 0.03 keV) while the reflected shock hard- 
ens the interior spectrum (/cTe = 0.12 to kT^ = 0.18 keV). The best spectral fits all required 
a depletion in oxygen. Attempts were made at nonequilibrium models but they did not im- 
prove the fit statistics. Furthermore, the fit values for the ionization parameter {uet, where 
He is electron density and t is the elapsed time since the gas was initially shocked) were all 
> 10^^ cm~^ s where values of > 3 x 10^^ cm~^ s signify equilibrium. The authors do not rule 
out the possibility of nonequilibrium but state that due to low spectral and spatial resolu- 
tion the low temperature equilibrium plasma is indistinguisha ble frorn a hig her temperature 
nonequilibrium plasma. Using 21 different extraction regions iLeahy I (120041 ) finds the same 
result; the spectra are best fit by an equilibrium MEKAL model with variable elemental 
abundances, and that there is no evidence for nonequilibrium. The author also states that 
the abundances are considerably depleted and vary on small spatial scales suggesting that 
the region is geometrically complex with multiple clouds and even more shocks. 



The final spectrum of note (IMiyata et al. 1120071 ) was taken using the Suzaku Observa- 
tory X-ray Imaging Spect r omete r (XIS) CCD camera and has the highest spectral resolution 
other than lVedder et al. I Jl986h . The authors clearly detected O VII at 562 ± 10 eV, O VIII 
at 653 ± 10 eV, C VI at 357 ± 10 eV and N VI at 425 ± 10 eV, but the 1/4 keV band 
emission is unresolved. Band ratio maps using strong emission features at different radii 
show that the outermost emission is dominated by the 1/4 keV emission and therefore a 
low temperature component. Furthermore, the O VIII/0 VII ratio increases inward showing 
that the ionization state is higher in the interior. The best fit models were two component 
nonequilibrium plasmas with variable abundances {vnei in Xspec). The lower temperature 
components of the fits varied from kTg = 0.10 — 0.15 keV and the high temperature com- 
ponent ranged from kTg = 0.18 — 0.34 keV. The temperature of both components increased 
towards the interior. Column densities ranged from = 3 — 6 x 10^° cm^^ and abundances 
were heavily depleted. 



The Cygnus Loop is clearly complex both spatially and spectrally. However, spectral 
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resolution is lagging spatial resolution. Chandra can image fine structures indicative of dif- 
ferent physical regions, but the nearly broadband spectra can only show trends. Determining 
plasma diagnostics from model fits to these low spectral quality data is uncertain. Higher 
spectral resolution is required not only to constrain the parameters of a model, but to test 
the assumed validity of the model. 

Currently there is no efficient, well developed technology that permits high resolution x- 
ray spectroscopy from large solid angle sources, making spectra of diffuse x-ray sources rare. 
In order to address this issue a soft X-ray spectrometer was designed for a rocket payload, 
the Cygnus X-ray Emission Spectroscopic Survey (CyXESS). Scientifically, the payload was 
designed to observe emission in the 1/4 keV bandpass. This is the least understood range of 
astrophysical soft X-ray energies. A diffuse high resolution spectrum has never been achieved 
even though there is a large amount of fiux in this band. The soft emission from the Cygnus 
Loop is contained within a narrow shell and is dominated by the early interactions of the 
initial blast wave with the surrounding cavity. Therefore, a global spectrum of the remnant 
will be dominated by the physics of this interaction. 



2. Sounding Rocket Instrument 



Typical X-ray telescopes employ the use of grazing incidence telescopes. However, these 
are expensive and heavy and thus unattractive for a rocket fiight. CyXESS utilizes a wire 
grid collimator to constrain the beam of light. The collimator is followed by the off-plane 
reflection grating array which disperses light onto the detectors. The payload is three meters 
long consisting of nearly a meter to create the converging beam while allowing the gratings 
to throw the light about two meters. A brief descr i ption o f the payload design follo ws, but 
a detailed description can be found in IMcEntaffer I (120071 ): IMcEntaffer et al. I ( 20081 ) . 



2.1. Wire Grid Collimator 

Wire grid collimators have wires that are spaced periodically in such a way that only 
light coming from a specified direction can pass through. If the grids have a spacing that 
decreases systematically, then it is possible to allow only light that is converging to a line 
to pass through, simulating the output of a lens. As light travels from front to back in the 
collimator it will encounter the same number of slits but they will be narrower and closer 
together, thus sculpting the converging beam. The wires create baffles between slits which 
vignette unwanted rays. If thin material is used, these wire spacers serve as knife edges 
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so that any light striking the metal will be near normal incidence and will be efficiently 
absorbed. Such a collimator does not function well for a point source, but for a diffuse 
target, radiation comes from all directions, and the beam is fully illuminated. A grating 
mounted in the exit beam diffiacts just as in a telescope beam. Thus the collimator alone 
provides the needed beam geometry. 



2.2. Off-plane Grating Array 

The off-plane mount at grazing incidence brings light onto the grating at a low graze 
angle, quasi-parallel to the direction of the grooves. The light is then diffiacte d through an 



arc, forming a cone , so that this mount is also known as conical diffiaction (jCash I Il991 



Catura et al. Ill988l ). The off-plane grating equation is 



smQ;-|-smp = - , (1) 

a sin 7 

where d is the spacing between grooves and 7 is the angle between the direction of the 
incoming ray and the direction of the groove at the point of impact. Light comes into the 
grating at an azimuthal angle of a along a cone with half-angle 7. It is then diffracted along 
the same cone of half-angle 7, but now with an azimuthal angle of f3. If 7 is kept small then 
the arc of diffraction stays close to the plane of the grating. However, at these low angles 
even large gratings have a small cross section to the incoming light. The solution is to place 
gratings in an array. 

The flight gratings have sinusoidal grooves with a density of 5670 ± 10 grooves/mm 
and an active area of 104 mm x 20 mm with the latter being the dimension along the 
groove. Using the limitations set by module packaging inside the payload as well as detector 
positioning on the bulkhead, the optimal graze angle is 4.4°. A 20 mm grating at 4.4° 
gives only a 1.53 mm effective length. Therefore, in order to maximize the throughput of the 
grating array, the grating substrates are extremely thin to minimize the amount of light they 
block. We used electroformed nickel with a figure of ~ A/2 and thickness of 0.127 ± 0.008 
mm. After accounting for the thickness, 67 gratings per module are required to catch the 
light in a 110 mm x 110 mm array. The array was held in a tension mount (~ 22 N per 
grating) to maintain the figure on the gratings and ensure survivability during launch. 
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2.3. GEM Detectors 

The two detectors on the payload are Gaseous Electron Multiphers (GEM). They were 
buih by Sensor Sciences, LLC. These innovative detectors use a gas filled chamber (75% 
Ar 25% CO2) segmented by perforated polyimide (Kapton(R)) film coated with a conductive 
layer on each side. The perforation holes provide a potential difference through which the 
electron cloud is accelerated resulting in gain. One of the most attractive features of these 
detectors is that they are made with very large formats, which is essential for this experiment 
due to the system's dispersion and line lengths. The entrance window is a 105 mm x 105 
mm polyimide window that is 3600-3900 A thick to maximize transmission while maintaining 
integrity. A 100 A carbon coat was added for conductivity and a grid bar and mesh support 
system is utilized with a transmission of 57.8%. The mesh and grid bars carry the negative 
high volts (HV) so that electrons are accelerated towards the anode, which is held at ground. 
Gain is determined by many factors including the voltage drop across the foils and gaps, 
high voltage supply stability, cleanhness of the GEM foils, gas pressure, etc. A gas flow 
system was incorporated to replenish the counter gas and maintain an operating pressure of 
14.5 psia. This system also counteracts the leak rate and compensates for micro tears in the 
window in order to improve gain stability. 

The 100 mm x 100 mm anode is a serpentine cross delay fine. The output of the resistive 
anode is analyzed by a custom electronics system. Signals are passed through an amplifier 
and then to an adder box which combines the data of the two detectors and passes it on to 
the timing-to-digital converter (TDG). The TDC returns a 12-bit word for X position and 
Y position and an 8-bit word for pulse height. The least significant bit of the pulse height 
determines which detector the current data word originates from. Finally, a stim pulse is 
sent from the TDG to the anode, which is then analyzed and sent back to the TDG. This 
gives reference data on both position and pulse height. This stim can be seen without HV 
on and thus provides a useful diagnostic. 

2.4. Expected Performance 

The payload was designed to obtain spectra from diffuse sources in the soft X-ray. 
Several design factors determine the accepted passband such as the length of the payload, 
the size of the detectors and the dispersion of the gratings. Optimizing these factors resulted 
in a passband of 44 A to 132 A in first order. The expected performance at these wavelengths 
can be summarized by the resolution and effective area of the payload. The resolution is 
ultimately determined by the full-width at half maximum (FWHM) of a spectral fine at 
the focus along with the dispersion of the system. The grating groove density of 5670 
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grooves/mm and throw distance of ~2 m give a dispersion of 0.89 A mm~^ in first order. 
Calibrations of the spectrum give hne widths that broaden shghtly from 1.7 mm to 2.2 mm 
as wavelength is increased. These characteristics give resolution of A/ A A ~ 25 — 70 in first 
order and ~ 25 — 85 in second order (22-66 A) as shown in figure [TJ 

The effective area of this spectrograph is determined by the collecting area, sky coverage, 
and throughput of the system. The collecting area of the telescope is defined by the size of 
the zero order image at the focal plane; 1.7 mm wide line over 10 cm of detector gives 1.7 
cm^. This small amount of collecting area is bolstered by the large amount of solid angle 
available to each point on the focal plane, 8.93 deg^. In terms of efficiency, the detector gas 
absorbs all X-rays, but the mechanical throughput is 57.8%, which is then further reduced 
by the transmission of the polyimide/carbon window. As for the gratings, the theoretical 
efficiencies are plotted as lines in figured Calibration data are shown as the points at 44.76 
A (carbon K-shell), with la Gaussian error bars, and agree well with theory. Taking all 
factors into account, the resulting effective area curves are shown in figure [31 

2.5. Flight 

The payload was launched from White Sands Missile Range at 02:00:00 UT, 2006 Novem- 
ber 21 (flight 36.224). The zenith angle of the target was ~ 31°. Usable data were recorded 
over 345 seconds of the flight. However, a breakdown event upon high voltage turn on ren- 
dered one detector useless while leaving the other detector extremely noisy for most of the 
flight. The GEM detectors exhibit noise in the form of hotspots which typically decay over 
time. Therefore, usable spectral data were only recorded over 65 seconds near the end of the 
flight. The pointing was dithered during the flight so that during the time when spectral 
data were collected the payload was pointed off center as depicted in figure HI 

3. Data Analysis 

Data were extracted from the section of the detector that was free of noise and contained 
spectral information. The resulting spectrum is given in figure O Pre-flight and post-flight 
spectral calibration data were compared with flight data to determine the wavelength scale. 
The counts at wavelengths longer than 50 A are residual emission from a large hotspot that 
dominated the detector for most of the flight. The two prominent features below 50 A are 
the detected spectral features. 

To fit the data a series of equilibrium spectral models of an optically thin plasma un- 
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der coUisional ionization e quilibrium were constructed using line lists from Raymond-Smith 

feavmond fc Smith Ill977h. MEK AL jMewe et al. Ill985l . Il98d : iKaastra lll992l:lArnaud fc RothenfluE 



1985 



Arnaud fc Raymond I Il992l : iLiedahl et al. I Il995f ) and APED fISmith et al. I l2001al Jbl) 

at temperatures ranging from kT^ = 0.034 — 0.272 keV (log(T) = 5.6 — 6.5). The spec- 
tra are constructed using Gaussians placed at the appropriate wavelengths (per Raymond- 
Smith, MEKAL or APED) with widths corresponding to the 1.84 mm (1.64 A) FWHM 
of the spectral line calibration data and amplitudes such that the integrated flux of each 
line scales according to the emissivities given in each line list. Each model spectrurn is ab - 
sorbed using photoelectric absorption cross sections from iMorrison fc McCammon I (119831 ) 



with Nh = 7 X 10 cm (held constant) and then conv olved w i th th e payload instrument 



response function (figure [3]). Cosmic abundances are per I Allen I (119731 ): He, 10.93; C, 8.52; 
N, 7.96; O, 8.82; Ne, 7.92; Na, 6.25; Mg, 7.42; Al, 6.39; Si, 7.52; S, 7.20; Ar, 6.80; Ca, 6.30; 
Fe, 7.60; Ni, 6.30 (in logarithmic units where log^Q Nh = 12.00). 

Given the low number of counts, a maximum likelihood analysis was performed to obtain 
the best fit. Assuming that the number of counts in each bin follows a Poisson distribution 
and that each bin is independent of the others results in the following likelihood function 



N 



i=l 



(2) 



The likelihood {£) was numerically maximized to provide the best fit and the most likely 
values for the parameters of the model (Zi) given the data (n^). The fit parameters include 
elemental abundances, which vary the amplitudes of corresponding lines, and a wavelength 
offset to the data, which adjusts the data wavelength scale within the errors of the detector 
calibration. Several priors were applied to constrain the fit parameters. First, abundances 
were assumed to always be > 0. Second, the wavelength offset (A shift) was always contained 
within the limits set by detector calibrations. Third, in all cases an attempt was made to 
vary as few abundances as possible while keeping abundances close to cosmic. Finally, fits 
were made with the constraint that line fitting only occurred at wavelengths less than ~50 
A while model lines longward of this were depleted until they were less than 2a above the 
long wavelength noise. Th e noise level was determined by a linear fit and the 2a upper limits 
were taken from iGehrels I (119861 ) for each bin given the fit. Only one model converged, a 
MEKAL plasma at kT^ = 0.14 (log(T) = 6.2). The fit to the spectral data was obtained by 
varying only the S abundance and the wavelength offset. The unintentional yet beneficial 
result that only S need be varied is not due to simplification but rather because S is the 
major contributor to the lines that need to fluctuate in order to obtain a reasonable fit. In 
order to reduce the model flux longward of 50 A to less than 2a above the noise. Si and N 
were depleted to 44% of cosmic. 
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The data and the model (sohd curve) are plotted together in figure [61 These data are 
identical to those shown in figure O but this time are plot ted with 0.9772 single-sided upper 
and lower limits (0.9544 confidence level) as calculated by lGehrels I ( 119861 ). which correspond 
to 2(j Gaussian statistics. The likelihood contours for the fit parameters are given in figure [71 
The shaded region encompasses 68% of the normalized likelihood and establishes the 84% 
(Icr Gaussian) marginalized confidence intervals for the individual parameters, 1. 55^0 63 for 
the S abundance and — 0.76^q;}7 A for the A shift. Finally, a closeup of the spectral data 
along with line identifications are shown in figure [HI The more prominent data line around 
44 A contains some first order Mg X and Si IX-Si XII, but most of the fiux is in the He-like 
triplet of O VII in second order. The other data line around 47-48 A is dominated by S IX 
and S X in first order. A summar y of the maj or lines is given in table [H The elemental 
abundances with respect to cosmic ( Allen 19731 ) are C, 1.0; N, 0.44 (max); O, 1.0; Ne, 1.0; 
Mg, 1.0; Si, 0.44 (max); S, 1.55+°;^°; Ar, 1.0; Ca, 1.0; Fe, 1.0; Ni, 1.0. Tabel □ summarizes 
the major lines and transitions. 



Discussion 



The results of the data analysis reveal a departure from cosmic abundances; S is enriched 
while Si and N are depleted. The enrichment of S can be explained by confusion due 
to multiple temperature com ponents. The single te mper ature fit at kT^, = .14 k eV is 



intermediate in comparison to iLevenson et al. I (Il999l) andlMivata fc Tsunemi I (120011 ) . but 
consistent with the low temperature component from lMiyata et al. I (120071 ). Therefore, there 
may be some contamination from a higher temperature component in the CyXESS data. 
In this passband, increasing the temperature increases the contribution from second order 
oxygen while decreasing the contributions of other lines. Therefore, a high temperature 
component will only add fiux to the line complex at 44 A, thus requiring additional fiux 
in the form of a S enhancement in the 47.5 A complex in order to maintain the spectral 
shape defined by the fit. Since there are only 2 features to fit, it is impossible to accurately 
discern the relative contribution from each plasma. More spectral resolution is required so 
that individual lines can be used to define the components. 

As for the second issue, even though some previous observations argue for an equilibrium 
plasma, especially in the case of the softest X-rays which occur in cloud shocks, nonequilib- 
rium conditions could still be important. If nonequilibrium is considered, then the ionization 
state of the gas should be decreased. In nonequilibrium the electron temperature is higher 
than what is expected from the ion state of the gas and the gas is underionized. Therefore, 
favoring lower ion species in the MEKAL model will test nonequilibrium conditions. Also, 
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within the lower ion species the ratio of higher energy hne flux to lower energy line flux should 
be increased relative to equilibrium since the electrons have more energy to excite the ions 
to higher levels than typical in equilibrium. In the case of Si, the lines that contribute the 
problematic long wavelength flux occur around 49-50 A and are Si IX, Si X and Si XI. De- 
creasing the influence of these lines will definitely have a desired effect, but will also require 
an increase of major Si VII and Si VIII lines at 52 A, thus reintroducing the problem. This 
occurs for other elements as well. Favoring lower ion species will not explain the depletions, 
and furthermore, there is no evidence for nonequilibrium conditions contributing to these 
data. 

Another explanation could be depletion into dust, especially in the case of Si depletion. 
If dust is a major constituent of the shocked ISM then refractory elements could be con- 
tained within the dust and depleted from the gas phase. However, the dust present in the 
cloud must be able to survive not only the initial blast wave, but also a reflection shock and 
sublimation over time in the hot gas (n ot to mention precursor w inds, cloud-cloud interac- 
tions, cosmic rays and photo desorption (Draine fc Salpeter 1979al lbl)). Given the conditions 
present in the Cygnus Loop, iDraine fc Salpeter I (1l979al ) show that thermal sputtering rates 
for dust grains in log(T) = 6.2 gas are ~0.001/im every 1000 years. Therefore, small grains 
(size <100A) will be quickly des troyed with the mass fraction returning to the gas. In addi- 
tion, iDraine fc Salpeter I (Il979bl ) have modelled dust sputtering as a function of blast wave 



velocity. Their results show that a blast wave with shock velocity > 300 km s~^ will sputter 
nearly all graphite, silicate, and iron dust grains up to a size of 0.1 fim in a cloud with 
density uh = 10 — 100 cm~^. Calculations for the in itial blast wave velocity for the C ygnus 
Loop vary from 330 km s~^ (ILevenson et al. II2002I ) to 400 km s~^ (IKu et al. Ill984l ). This 
also suggests that if a significant fraction of refractory elements are depleted into dust in 
the dense clouds consituting the Cygnus Loop cavity walls, then the size distribution of dust 
particles favors large grains. 



Infrared (IR) emission from the Cygnus Loop has been observed and lArendt et al 



(I1992I ) show that it can be explained by dust emission. Using IRAS observations at 12 fim, 
25 /im, 60 /xm, and 100 /xm, they find that there are two infrared (IR) components, one that 
correlates well with X-ray emission and another that correlates with the optical emission. 
A lack of observed emission in the 12 /xm and 25 fim bands suggests an underabundance of 
small grains. This is supported by their models of the X-ray/IR correlated gas, which favor 
a minimum grain size of ~150 A consistent with our estimates. However, the models used 
to fit the broadband observations assume all emission is due to thermal dust emission. The 
authors state that line emission could theoretically contribute a significant fraction to the 
observed IR emission. Therefore, constraining the dust fraction is impossible without higher 
resolution IR spectra, which should become available with upcoming Spitzer observations. 
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Depletion into dust cannot be ruled out and is supported by our data, especially in the 
case of silicate grains. We see no evidence for graphite grains because there are no important 
C lines in the wavelength range of our spectrum. The N and Si depletions necessary for our 
fit are due to a line complex of second order N VII and first order Si X and Si XI. Therefore, 
given this confusion only an upper limit to the Si abundance can be applied suggesting a 
large depletion of at least 56% of the Si into grains. 

5. Summary 

The CyXESS payload was designed to observe the soft X-ray fiux of the Cygnus Loop 
supernova remnant. The design consisted of a wire grid collimator that focused the light 
onto an array of gratings in the off-plane mount which ultimately dispersed the spectrum 
onto large format GEM detectors. The payload was launched on November 20th, 2006 from 
White Sands Missile Range and collected 345 seconds of data. Data reduction decreased 
the amount of usable data to 65 seconds during which the instrument detected flux between 
43-49.5 A (250-288 eV) in two prominent features. The flrst feature near 45 A is dominated 
by the He-like triplet of O VII in second order with contributions from Mg X and Si IX- 
Si XII in first order, while the second feature near 47.5 A is first order S IX and S X. Fits to 
the spectra give an equilibrium plasma at kT^ = 0.14 keV (log(T) = 6.2) and near cosmic 
abundances for most elements. Even though the most likely fit to the CyXESS data contains 
only one temperature, the lack of spectral range and resolution do not allow determination 
of multiple components. An observed depletion in Si supports the presence of silicate grains 
but higher X-ray and IR resolution are necessary to accurately constrain dust models. 

Our data were constrained to a small portion of the soft X-rays, but show a wealth 
of lines present. In this band we find no evidence to support nonequilibrium conditions. 
Therefore, our results are consistent with previous observations and show that the soft X- 
ray spectrum of the Cygnus Loop, which is dominated by interactions between the initial 
blast wave with the walls of a precursor formed cavity, can be described by an equilibrium 
plasma. 

The authors would hke to acknowledge NASA grants NNG04WC02G and NGT5-50397 
for support of this work. 
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Fig. 1. — Spectral resolution as a function of wavelength for first order (solid line) and second 
order (dashed line). 
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Fig. 2. — Theoretical grating efficiency curves for first order (dashed), second order (dot- 
dashed) and sum of orders (sohd). The data points are from pre-fiight cahbrations using 
carbon K-shell emission. 
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Fig. 3. — Effective area curves for first order (dastied), second order (dot-dastied) and sum 
of orders (solid) . 
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Fig. 4. — The octagonal FOV of CyXESS is sho wn over a ROSAT brightness map of the 
0.25 keV band taken from iLevenson et al. I (119991 ) . This pointing is ~ 55' to the west of the 
center position and is where the telescope was pointed when the spectral data were observed. 
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Fig. 5. — Two spectral lines stand out on the left side. The rest of the counts are residual 
broad oval emission that dominates the noise. 
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Fig. 6. — Data are plotted as points with 2(7 errorbars. The log(T') = 6.2 MEKAL equihb- 
rium plasma model is plotted as the solid line. 
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Fig. 7. — Likelihood contours for the fit parameters are plotted in steps of 10% of the peak 
value. The gray area encompasses the 68% confidence region. 
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Fig. 8. — Line identifications for the spectral data. At the top, the ion names point to 
vertical lines that depict their relative contributions. Second order lines are called out above 
the first order hues. Spectral data are plotted as the histogram. The MEKAL plasma fit 
is plotted as diamonds with first order model contributions as asterisks and second order 
contributions as diamonds. 
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Table 1: Major spectral lines (letected ])\ C}-XESS. 



Ion 


Wavelength (A) 


Transition 


Relative strength 


VII 


21.600 


ls2 - Is 2p [R] 


1.000 


VII 


21.800 


Is^ - Is 2p [I] 


0.216 


VII 


22.100 


ls2 - Is 2s [F] 


0.852 


Si XI 


43.740 


2s2 - 2s 3p 


0.182 


Mg X 


44.050 


2s - 4p 


0.116 


Si XII 


44.165 


2p- 3d 


0.096 


SIX 


47.500 


2p4 - 2p3 3d 


0.444 


sx 


47.654 


2p3 - 2p2 3s 


0.378 


sx 


47.793 


2p3 - 2p2 3s 


0.186 



